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■ Abstract: We define stabilizability of an infinite volume height configuration and of a 
, probability measure on height configurations. We show that for high enough densities, 

— , ■ a probability measure cannot be stabilized. We also show that in some sense the 

thermodynamic limit of the uniform measures on the recurrent configurations of the 
abelian sandpile model (ASM) is a maximal element of the set of stabilizable measures. 
In that sense the self-organized critical behavior of the ASM can be understood in terms 

■ of an ordinary transition between stabilizable and non-stabilizable. 
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>$ : 1 Introduction 

a ■ 

Self-organized criticality (SOC) is a concept introduced in [2j to model power-law 
behavior of avalanche sizes in various natural phenomena such as sand and rice piles, 
^ ■ forest fires, etc. The conceptual point of view in [Hj is that this kind of criticality 
is not tuned by parameters such as temperature or magnetic field, as is the case in 
critical systems of equilibrium statistical mechanics. This point of view has been 
questioned by several people, see e.g. 0111, where it is argued that the choice of the 
models exhibiting SOC involves an implicit tuning of parameters, and hence SOC is 
an (interesting) example of ordinary criticality. In the case of the abelian sandpile 
model, e.g. one can say that the choice of the toppling matrix (which governs the 
dynamics) having mass equal to zero is a fine tuning. Indeed in the massive or 
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dissipative case (where in the bulk grains are lost upon toppling) the avalanche sizes 
exhibit exponential decay, so in that case there is no criticality. 

Similarly, in [llj . the authors investigate the relation between the critical density 
of some parametric model of random walkers with that of the abelian sandpile model, 
and prove in d = 1 that ASM density corresponds exactly to the transition point in 
the random walkers model. They further conjecture that this is also true in d > 2. In 
this paper we want to continue the relation between an ordinary critical phenomenon 
and the SOC-state of the abelian sandpile model. This is done through the notion of 
"stabilizability" . A height configuration is called stabilizable if upon stabilizing it in 
larger and larger volumes, the number of topplings at a fixed site does not diverge. 
This implies that we can "redistribute" the mass in infinite volume such that after 
the redistribution, all sites have a height between 1 and 2d. Similarly a probability 
measure v on height configurations is called stabilizable if it concentrates on the set of 
stabilizable configurations. The conjecture in [TT], inspired by [S] is that there exists 
p c > such that (modulo some restrictions on the measure v) if the v expected height 
p < p c , then v is stabilizable, if p > 2d it is not stabilizable, and for any p G (p c , 2d) 
there exist measures v with expected height p which are not stabilizable. Moreover, 
p c is exactly the expected height in the stationary state of the abelian sandpile model, 
in the thermodynamic limit. 

The aim of this paper is to prove the last two items of this conjecture, and to 
give some more insight in the regime p G (p c , 2d). Our paper is organized as follows: 
in section 2 we define the notion of stabilizability. In section 3 we precisely state 
the main conjecture of JTJj and prove item 3 of it. In section 4 we prove item 2 of 
the conjecture and disprove by a counterexample item 1. We further show that in 
some sense the infinite volume limit of the stationary measure of the abelian sandpile 
model is "maximal stabilizable". Finally in section 5 we introduce the concept of 
meta-stabilizability, and give a class of examples of probability measures on height 
configurations having this property. 

2 Basic definitions 

A height configuration is defined as a map rj : 7/ — > {1,2,...}, where by convention 
the minimal height is chosen to be 1. The set of all height configurations is denoted 
by TC. The set of probability measures on the Borel sigma- field of TC is denoted 
V(TC). A configuration is called stable if for all x G Z d , rj(x) < 2d. The set of stable 
configurations is denoted by Q. Similarly, for V C Z d , Qy denotes the set of stable 
configurations 77 : V — > {1, . . . 2d}. 

For V G 7j d the abelian sandpile toppling matrix is defined as 



(Ay) 




(2.1) 



Its inverse is denoted by 



G v (x,y) = (Ay 1 ) 



(2.2) 
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The probabilistic interpretation of Gy is: 

Gy(x,y) = — Z 1 ,. (number of visits at site y) (2.3) 

where E]f denotes expectation in the simple random walk started at x, and killed 
upon exiting V. 

Definition 2.4. A height configuration r\ G Ti. is called stabilizable if for any sequence 
of volumes V n f Z d ; there exist my n G N v " such that 

Wn - A Vn m Vn = £y n G VL Vn (2.5) 

and for any x G 7L d my n (x) — » m(x) as n — > oo. 

The set of all stabilizable configurations is denoted by S. It follows immediately 
from the definition that for r\ G S, and m = lim n ^ 00 , 

ij-Am = (G(] (2.6) 

where A XjV = 2d5 xy — lu-^i is the infinite volume toppling matrix. 

Definition 2.7. A probability measure v on height configurations is called stabilizable 
»/!/(£) = 1. 

It is clear that the set of stabilizable configurations is a translation invariant subset 
of 7i. Therefore any stationary and ergodic probability measure \i on TC satisfies 
/x(5)€ {0,1}. 

Remark that in finite volume V C Z d , the abelian sandpile model is "well-defined". 
This means that for any height configuration rj G N y , the equation 

7] V - A v m v = £ y (2.8) 

with unknowns the couple (my, £y) has at least one solution, namely for x G V, 
mv(x) equals the number of topplings at x needed to stabilize r\ in V (see e.g. H). 
Notice that the couple (my, £y) is not unique, but if my is such that £y is stable and 
()2.8|) holds, then m v (x) is at least the number of topplings at x needed to stabilize 
rjy. So the vector collecting the number of topplings (needed to stabilize r/y) is the 
minimal solution my of the equation ()2.8j) . We will always choose this solution in the 
sequel. 

On Ti we have the natural pointwise ordering rj < £ if Vx G Z , r/(x) < ^(x). 
A function / : 7^ — > IR. is called monotone if rj < £ implies /(^) < /(£). Two 
probability measures fi,v on H, are ordered // < i/ if for any monotone function, 
/ /^A* — / fdv- This is equivalent to the existence of a coupling P of fi and v such 

that P({(tj, £}) = !■ 

By abelianness, my(x) is non-decreasing in V. Therefore a configuration 77 is not 
stabilizable if and only if there exists x G Z d such that my(x) j 00. By abelianness, 
the my are also monotone functions (in the sense just mentioned) of the configuration 
V- 

Therefore we have the following immediate properties of the set of stabilizable 
configurations 
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Proposition 2.9. a) S is a translation invariant measurable set. 

b) If i] G S and £ < r\, then £ G S. 

c) If fi is a stabilizable probability measure, and v < p,, then v is a stabilizable 
probability measure. 

We then define the following "critical densities": 

Lemma 2.10. Define 

pi = inf{p > 1 : 3v G V{H) with v(r)(0)) = p, and v is not stabilizable} 

p~ = sup{p > 1 : Wu G V(H) with z/(?7(0)) = p, v is stabilizable} (2-11) 

Then p+ = p~ . 

Proof. It suffices to see that the set 

S = {p > 1 : such that \/v G V{H) with v(r)(0)) = p, v is stabilizable} (2.12) 

is an interval. Suppose that p G S and p' < p. Consider a measure v' G V(H) 
with v'(r)(0)) = p' . Then there exists a measure v G P('H) such that v(r)(0)) = p 
and v > v' . Since v is stabilizable, by the monotonicity property 12.91 item 3, v' is 
stabilizable. □ 

We now introduce the "critical state" of the sandpile model, and its thermody- 
namic limit. Define a configuration allowed in a volume V if for any subset W C V, 
the inequality 

v(x)<\{yeW,\y-x\ = l}\ (2.13) 

is violated for at least one x G W . The set of allowed configurations in volume V is 
denoted by IZy. It is well known that the set of recurrent configurations of the abelian 
sandpile model in finite volume V coincides with the set of allowed configurations TZv, 
and the stationary measure is the uniform measure on IZy- (see e.g. jH] or the basic 
reference jl]). We denote this measure by p,y. Recently, it has been proved in 12] 
that the weak limit p = limyf Z d \iy exists and defines a measure on infinite volume 
height configurations. Moreover, its support 1Z is the set of those configurations rj 
such that all restrictions to finite volumes V have the property r\y G IZy. We will call 
this measure \x the uniform measure on recurrent configurations (UMRC). We will 
always use the symbol fi for the UMRC and denote its mean height by p c = /x(r/(0)). 
In [2] it is proved that \x is translation invariant, while in [Hj it is proved that p, is 
tail-trivial so in particular ergodic under spatial translations. 
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3 Main conjecture and results 

In the rest of the paper we will prove point 2 and 3 of the following conjecture 
appearing in [TT], cf. also [5], and we will also give some additional new results and 
examples. 
Conjecture: 

Let v be a stationary and ergodic probability measure on N z . Put p = 1/(77(0)). 

1. For p < p c , v is stabilizable 

2. For p c < p < 2d there exist v which are not stabilizable 

3. For p > 2d, v is not stabilizable 

The following theorem settles point 3 of the conjecture. 

Theorem 3.1. Suppose that rj has a distribution v such that 1/(77(0)) = p > 2d. Then 
7] is almost surely not stabilizable. 

Proof. Let Gy be the Green function introduced in (j2.1|) . ()2.2|) . Simple random walk 
killed upon exiting V will be denoted by {X n , n e N}, and corresponding expectation 
by Kq . Finally, let Ty denote the lifetime of this walk. The infinite volume random 
walk expectation is denoted by E . Of course, E and Ejf-expectation of events before 
Ty coincide. 

Suppose that 77 drawn from v is stabilizable. Then we have 



and my(0) j m(0) < 00 as V j 7L d . Since the random field 77 is stationary and ergodic, 
we have 



Po x v almost surely, where Po denotes the path-space measure of the simple random 
walk starting at 0. For £ we cannot conclude such a strong statement but we have, 
by stability 



m v {0) 



]TGV(0,x)(77(x)-£(x)) 




(3.2) 




limsup — \ £(X n ) < 2d 
v\z d T v n=1 

Therefore since p > 2d + S for some 5 > 0, 





(3.3) 
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This implies, using Fatou's lemma, and the fact r v — > oo as V | Z d , that for any 
A > 

/ TV 

2dliminfm y (0) = liminfE V(n(X n ) - £(X n ) 

V\1 d V\TL d > 



v n=0 



> lim inf E Ah Tv>A) 



> 



> (3.4) 
Since A > is arbitrary, we arrive at a contradiction. 

□ 



4 Adding to the stationary measure 

In this section we settle point 2 of the conjecture. 

The UMRC fx is obtained as a limit of finite volume stationary measure fiy. These 
fiy are in turn obtained by running the finite volume addition and relaxation process 
for a long time. Therefore, one can believe that fx is "on the edge" of stabilizability. 
More precisely if one could still "add mass" to fx, then fx would not be stationary. 

However, it is not true that fx is a maximal stabilizable measure in the sense of the 
FKG ordering of measures. Indeed, one can create the following translation invariant 
v: pick a configuration according to fx and flip all the height ones to height four. 
This measure is strictly dominating fx in FKG sense, but it concentrates on stable 
configurations. In the last section of this paper we will show that such "artificially 
stable" measures are in some sense "metastable" . 

The idea of formalizing the maximality of fx is that "one cannot add mass to fx" . 
For fx a probability measure on Tt, and v a probability measure on N zd , we denote 
by fx © v the distribution of r\ + £ where r\ is distributed according to fx and £ is 
independent of rj and distributed according to v. 

Definition 4.1. A probability measure fx on 7i is called maximal stabilizable if for 
any ergodic translation invariant v with v(r)(0)) > 0, fx © v is not stabilizable. 

In order to state our main result of this section, we need some more conditions 
on the UMRC fx. For a configuration drawn from fx, define the addition operator 
a x,v by a x y{rj) = (a x yr] V )r]y c . In words this means that application of a x y to a 
configuration in V, te configuration changes as if we added at x and stabilized in V, 
while outside V the configuration remains unaltered. 

We say that the infinite volume addition operator a x = lim v ^ Z d(a x y) is well- 
defined w.r.t. the UMRC if for fx almost every rj, the limit lim^^d a x y(rj) exists (in 
the product topology). We now can state our conditions 
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Definition 4.2. The UMRC is called canonical if 

1. The infinite volume addition operators a x are well-defined w.r.t. the UMRC for 
any x G Z* d . 

2. The UMRC is stationary w.r.t. the action of a x , i.e., ifrj is distributed according 
to the UMRC, then so is a x rj. 

In jHj we prove that these conditions are satisfied on Z d , d > 5. The restriction 
d > 5 is however of a technical nature, and we strongly believe that these conditions 
are satisfied as soon as \i exists. If the UMRC fi is canonical, then one can easily see 
that finite products of addition operators are well-defined fi a.s. and leave fi invariant. 
See |H] for a complete proof. 

Our main result in this section is the following. 
Theorem 4.3. If the UMRC is canonical, then it is maximal stabilizable. 

Proof. We have to prove that fx@v is not stabilizable for any v stationary and ergodic 
such that 1/(77(0)) > 0. A configuration drawn from \i © v is of the form 77 + a, where 
77 is distributed according to fi and a independently according to v. 

Suppose 77 + a can be stabilized, then we can write 

r] V + a v - A v m v = £ y (4.4) 

with m\ | m\ d as V | Z d . We define m 2y G N zd by 

77 + a° v - Am 2y = i 2y (4.5) 

where a v : Z d — >• N is defined a v (x) = a^x)!^^- in words this means that we add 
according to a only in the finite volume V but we stabilize in infinite volume. The 
fact that m 2,v is finite follows from the fact that the addition operators a x and finite 
products of these are well-defined in infinite volume on /1 almost every configuration. 
Since for W D V 

a° v < a° w (4.6) 

and Tfiy does not diverge, it is clear that m 2,v is well-defined, by approximating 
the equation ()4.5|) in growing volumes. Moreover, for A C Z d fixed, it is also clear 
that (m 2 ' y )A and (my) a will coincide for V D Vq large enough. Otherwise, the 
stabilization of 7/y + ay would require additional topplings in A for infinitely many 
V's, which clearly contradicts that m v converges (and hence remains bounded). But 
then we have that for V large enough, (£ v )a and ^ v coincide. For any V, the 
distribution of £ 2 ' y is //, because \i is stationary under the infinite volume addition 
operators. Therefore, we conclude that the limit limy £ v = limy £y is distributed 
according to /1. Hence, passing to the limit V f Z d in (|4.5|) we obtain 

77 + a - Am = f (4.7) 
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where r\ and £ have the same distribution /i, and where m G N zd . Let {X n ,n G N} 
be simple random walk starting at the origin. Then for any function / : Z d — > R, 

n— 1 

M n = /(X n ) - /(X ) - ^ ^(-A)/(X i ) (4.8) 

i=i 

is a mean zero martingale. Applying this identity with f(x) = m x , using ()4.7|) gives 
that 

n-1 

M n = m(X n ) - m(X ) + ^(r/(X 4 ) + a(X t ) - £(X t )) (4.9) 

i=i 

is a mean zero martingale (w.r.t. the filtration T n = a(X r : < n), so 77 and £ are 
/ixed here) . 

. This gives, upon taking expectation over the random walk 

- (E (m(X n )) - m(0)) = -E f Y> (X*) - ^(X,) - «(X,)) ] (4.10) 



i=i 



Using now that m(0) < oo by assumption, the ergodicity of \i and u, and the fact 
that both £ and 77 have distribution /x, we obtain from (|4.10j) upon taking the limit 
n — >• 00 that 



—a 



< lim i (E (m(X n )) - m(0)) = lim E (- V(£(X,) - r/(X,) - a(X l )) ] 

\ i=l / 

(4.11) 

which is a contradiction. □ 
In dimension d = 1 the situation is simpler, see fT] for the proof. 

Theorem 4.12. A stationary and ergodic measure v on N z such that i/(j](0)) < 2 is 
stabilizable. If on the contrary i/(j](Q)) > 2, then v is not stabilizable. 

Remark 4.13. For is(rj(0)) = 2 one can have both stabilizability and non-stabilizability: 
e.g. the configuration 313131313 . . . and its shift 13131313 . . . are not stabilizable. 



5 Constructive example 

In this section we present an explicit example of an addition which leads to infinitely 
many topplings at the origin in the limit V | It settles point 2 of the conjecture, 
even in the case when the UMRC is not canonical (in particular for d < 4), and 
shows that point 1 is not true in the generality of stationary and ergodic probability 
measures on height configurations. 

Before presenting the example, let us recall the following fact about the abelian 
sandpile model, see e.g. jTUj. Consider the abelian sandpile model in a finite volume 
V C 7j d . Start from a recurrent stable height configuration r\ : V — > {1, . . . , 2d}. Add 
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on each boundary site x G dV as many grains as there are "lacking neighbors", i.e., 
the number of grains added at x equals Ay(x) = |{?/ 6 Z d : \y — x\ — 1, y ^ V}|. Then, 
upon stabilization, each site will topple exactly once and the recurrent configuration 
7] will remain unaltered. If V is a rectangle, and a stable height configuration in 
V is recurrent, then this "special addition" consists of adding two grains to the 
corner sites and one grain to the other boundary sites. If rj : V — > {1, . . . , 2d} is 
recurrent in V, and W C V, then the restriction r]w is recurrent in W. Therefore, 
by abelianness, if we add to each boundary site x of W at least a number grains 
= {y £ Z d : \y — x\ — l,y ^ W}, then, upon stabilization, each site in W will 
topple at least once. 

We can now present our example in the case d = 2; the generalization to d ^ 2 is 
obvious. Let u,u' be independent and distributed according to a Bernoulli measure 
P p on {0, 1} Z , with F p (uj(x) = 1) = p. Consider the following two dimensional 
random field ((x, y) = u(x) + w'(y). This is what we are going to add to a recurrent 
configuration. In words, if for x G Z, u(x) = 1, then we add one grain to each lattice 
site of the vertical line {(x, y), y G Z}, and if u'(y) = 1 then we add one grain to each 
lattice site of the horizontal line {(x, y) : x G Z}. If we add according to (, then there 
are almost surely infinitely many rectangle Ri, . . . R n , . . . surrounding the origin with 
corner sites where we add two grains and other boundary sites where we add at least 
one grain. 

If we add such a configuration ( to any recurrent configuration rj, then we have 
that the number of topplings at the origin in the finite volume V is at least the number 
of rectangles Ri that are inside V. Indeed, upon addition according to ( on such a 
rectangle, every site inside the rectangle will topple at least once. 

Therefore the distribution p p of rj + ( where i] is drawn from the UMRC /i, is not 
stabilizable. Since we can choose p arbitrary close to zero, any density p G (p c , p c + 2) 
can be attained by fi p . 

To show that we can actually get below p c , remember that the fact that the number 
of topplings inside V is at least the number of rectangles Ri does only depend on the 
fact that the configuration to which we add is recurrent, and not on the fact that the 
configuration is chosen from a particular distribution. 

Therefore, consider a translation invariant probability measure // concentrating 
on a subset TV of 1Z. Then with the same reasoning, the distribution fi' p of rj + ( where 
7] is drawn from the /i', is not stabilizable. Consider therefore // to be a weak limit 
point of the uniform measures on minimal recurrent configurations, where "minima" 
is in the sense of the pointwise ordering of configurations. Then the distribution p' p 
has its expectation yU p (^(0)) arbitrary close to //(//(O)). This number is strictly less 
than the critical density p c . Indeed, it is at most 3, because the all three configuration 
is recurrent, and p c > 3, see e.g., [T2] . 

This shows that point 1 of the conjecture cannot hold in that generality. 

Combining our results so far with proposition 2.4 from [TI], we conclude 

Theorem 5.1. Let p+ = p~ be as in lemma \Q.l(A Then 

p~^ = inf{i/(r](0)) : v is translation invariant and v{TV) = 1}. (5-2) 
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For any dimension d, p+ < p c , and for d = 2, p+ < 3 < p c 

Remark 5.3. We believe that the strict inequality p+ < p c holds in any dimension 
d > 1, but this is an open problem as far as we know. In d = 1 we have p+ = 2 = p c , 
but this is an exceptional case where almost all recurrent configurations are minimal 
recurrent. 

6 Other notions of stabilizability 

6.1 Stabilization in infinite volume 

Definition 6.1. A configuration rj G N zd is called weakly stabilizable if there exists 
meN zd and £ G O = {1, . . . , 2d} zd such that 

r? - Am = £ (6.2) 

It is clear that if rj is stabilizable, then it is weakly stabilizable and we can choose 
m = liniy| Z d my. However it is not clear whether there exist unstable configurations 
which can be stabilized directly in infinite volume but which satisfy limyrad my(0) = 
oo, i.e., the infinite volume toppling numbers are not obtained as the limit of toppling 
numbers in larger and larger volumes. In the following proposition we prove that a 
measure v with v{r\{xj) > 2d for all x G Z d cannot be weakly stabilized. This means 
in words that mass cannot be "swept away" to infinity. 

The following example shows that the opposite, importing mass from infinity, is 
not impossible. Consider / : Z 2 — ► N: 

f(x,y) = x 2 + y 2 

then Af = —4 and hence, for example, 

6 = 2 — A f (6.3) 

where 6 (resp. 2) denotes the configuration with height 6 (resp. 2) at every site. 
However, 

lim A^(6 - 2) = oo 

viz d 

so this "infinite volume toppling" cannot be obtained as a limit of finite volume top- 
plings. Notice however that "toppling" according to / is not "legal" in the following 
sense: we cannot find an order of topplings such that, performed in this order, only 
unstable sites topple and at the end every site has toppled x 2 + y 2 times. The point 
of ()6.3j) is that the equality 

T] = £ — Am (6-4) 

for some m G N zd does not imply that the densities of r\ and £ are equal. However as 
we will see later, the equality in ()6.4j) does imply that the density of r\ is larger or 
equal than that of £. 

In the following proposition we show point 3 of the conjecture for weak stabiliz- 
ability. 
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Proposition 6.5. Let v be stationary ergodic such that 1/(77(0)) > 2d. Then v is not 
weakly stabilizable. 

Proof. Suppose that there exist m G N z<i such that 

77 — Am = £ (6.6) 

with £ stable and 77 a sample from v. Let X n be the position of simple random walk 
starting at the origin at time n. ^From ()6.6j) it follows that 

- n— 1 

m(X„) - m(0) - - $>(X fc ) - 7/(X fc )) (6.7) 

k=l 

is a mean-zero martingale. Therefore taking expectations w.r.t. the random walk 

1 1 / 1 n_1 \ 

- (E (m(X n )) - m(0))) = -E - - rj(X h )) (6.8) 
n n \ 2d 4 ' 

\ fc=i 

By stability of £ 

- ( - itf**)) ) < 2d - \ 5>pQ) (6.9) 



n—l \ ., n— 1 



wfc=l / fc=l 



Therefore, using dominated convergence and z/(r/(0)) = 2d + 5, 
< liminf - (E (m(X n )) - m(0)) 

1 / 1 

< hmsup -E — ^(£(X fc ) - rj(X k )) 
n— too n \ zd 



< l-Eof-L lim -J2v(Xk) J <S (6.10) 
\ fe=i / 



which is a contradiction. 



□ 



6.2 Activated random walkers system and stabilizability at 
low density 

Dickman proposes in [Sj the following mechanism of stabilization. Consider a config- 
uration 77 G TC. To each 2 G Z d is associated a Poisson process N£, for x ^ y these 
processes are independent. On the event times of Nf a site topples if it is unstable 
(the random walkers are "activated"), otherwise nothing happens. This means that 
after time t the configuration 77 evolves towards r\ t according to a Markov process with 
generator 

Lf( V ) = £ l v{x)>2d (f( V - A Xj .) - f(rj)) 

x&L d 
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One says now that the configuration is stabilizable by this process if for any x G Z d 
the value r) t (x) jumps only a finite number of times. One can write the configuration 

r)t as 

Vt = Vo~ An* (6.11) 
where n* is the vector collecting at each site the number of topplings at x in [0,t]. 
If 7] is distributed according to a translation invariant probability measure v on N zd 
then so is under the joint measure v x P where P is the distribution of the Poisson 
processes iV*, x G Z d . Moreover [x] < Nf by definition and hence 

E( Vt (x)) = E( Vo (x)) 

i.e., this process conserves the density. 

Lemma 6.12. A configuration is stabilizable by the process with generator L if and 
only if it is stabilizable (in the sense of definition \2. 



Proof. Suppose r\ is stabilizable by the process with generator L. Consider then the 
generator 

Lvfiv) = E^)>m W 7 ? " A £) - fW) 

x&V 

corresponding to toppling inside V only, according to the finite volume toppling ma- 
trix, on the event times of the Poisson process iV*. Call n v (x, rj) the number of 
updates of x in [0,t]. It is easy to see that n v (x, rf) f n\.{x) as V | Z d , where n* is 
defined above. By assumption, for any x there exists t^(rj) such that for any t > t^(rj) 

n v (x,f]) = riy ^ v \x,r]), and moreover t^(j]) j t x (r]) for which n*(x) = n^ v (x) for any 
t > t x (i]). Therefore in definition 12 .71 we can identify 

my (x) = Vby ^ (x, 7]) 

and 

m(x) = n^^x) 

Suppose that 77 is stabilizable in the sense of definition 12.71 Then, clearly, in finite 
volume we have the equality 



Since my fmwe have 



and hence 



my{x) = n-y (x, 7]) 



sup n-y [x, r]) < oo 

vci d 



tx{v) = sup tl(ri) 

vcz d 



is finite P almost surely. Now pick t > t x (r]). Then 



nAx^n) = lim ny(x,ri) = lim ny^(x,ri) (6.13) 
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where in the second step we used that the processes with generator L v converge to 
the process with generator L weakly on path space. Indeed for any local function 

lim L v (f) = L(f) 

So the convergence of the processes follows from the Trotter-Kurtz theorem. The 
right hand side of (j6.13|) does not depend on t anymore. Hence 7] is stabilizable by 
the process with generator L. □ 

In [TT] the authors prove that there exists p' c > such that if z/(^(0)) < p' c , then v 
is stabilizable by the process with generator L. p' c is the density of "minimal recurrent 
configurations". The following theorem is then an immediate consequence. 

Theorem 6.14. There exists p' c > such that ifv is any stationary ergodic probability 
measure on 7i with ^(?](0)) < p' c , then v is stabilizable. 

Proof. Combine lemma IB . 1 21 with proposition 2.4 from jTTj. □ 



7 Metastable measures 

Suppose that v > p, and v(r)(0)) > p(r](0)), i.e., v has a strictly higher density than 
p, and stochastically dominates p. In that situation v can still concentrate on stable 
configurations, and hence be stabilizable. One feels however that such a measure is 
"on the brink of non-stabilizability" . This is formalized in the following definition. 

Definition 7.1. A measure v is called metastable if it is stabilizable and if u © 5q is 
not stabilizable with non zero probability, i.e., v © S (S) < 1. 

In words this means that upon stabilizing v © 5q in volume V, with positive prob- 
ability the number of topplings my(0) diverges as V t The simplest example of 
a metastable measure is the measure concentrating on the maximal stable configura- 
tion v = #2d- The following theorem shows that there are other non-trivial metastable 
measures. 

Theorem 7.2. Suppose that v is a stationary and ergodic probability measure on Q, 
concentrating on the set of recurrent configurations TZ. Define I v {x) = \( x )=2d and 
call v the distribution of I v . Suppose that v dominates a bernoulli measure F p with 
p sufficiently close to one such that the 1 ; s percolate and the zeros do not percolate. 
Then v is metastable. 

Before giving the proof, we need the concept of "wave", see e.g. [TU]. Consider 
a stable height configuration in a finite volume V. Add to x and topple x once (if 
necessary) and perform all the necessary topplings except a new toppling at x. This 
is called the first wave. The support of the wave is the set of sites that toppled (it 
is easy to see that during a wave all sites topple at most once). If x is still unstable, 
then iterate the same procedure; this gives the second wave etc. For the proof of 
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theorem 17. 2[ we need that the support of a wave is simply connected, i.e., contains 
no holes. For the sake of completeness, we give a precise definition of this and prove 
this property of the wave, already formulated in [TO] . 

Definition 7.3. A subset V C 7L d is called simply connected if the set obtained by 
"filling the squares", i.e., the set 

U, ey (* + [-l/2,l/2] d ) 

is simply connected. 

Lemma 7.4. Let r\ G IZy be a recurrent configuration. Then the support of any wave 
is simply connected. 

Proof. Suppose the support of the first wave contains a maximal (in the sense of 
inclusion) hole ^ H C V. By hypothesis, all the neighbors of H in V which do 
not belong to H have toppled once. The toppling of the outer boundary gives an 
addition to the inner boundary equal to Xh{x) at site x, where Xh(x) is the number 
of neighbors of x in V, not belonging to H. But addition of this to t]h leads to one 
toppling at each site x G H, by recurrence of T)h- This is a contradiction because we 
supposed that H ^ 0, and H is not contained in the support of the wave. After the 
first wave the restriction of the configuration to the volume V \ x is recurrent i.e., 
recurrent in that volume V \x. Suppose that the second wave contains a hole, then 
this hole has to be a subset ofV\x (because x is contained in the wave by definition), 
and arguing as before, one sees that the subconfiguration r] H cannot be recurrent. □ 

We can now give the proof of theorem ([7.2)1 . 

Proof. The idea of the proof is the following. Suppose we have a "sea" of height 2d 
and "islands" of other heights, and such that the configuration is recurrent. Suppose 
the origin belongs to the sea, and we add a grain at the origin. The first wave must 
be a simply connected subset of Z d because the configuration is recurrent. It is clear 
that the "sea" of height 2c? is part of the wave, and therefore every site is contained 
in the wave (because if an island is not contained then the wave would not be simply 
connected). So in the first wave every site topples exactly once, but this implies that 
the resulting configuration is exactly the same. Hence we have infinitely many waves. 

Let us now formalize this. For a given configuration T) (distributed according to 
v) a volume V C Z d is called "a lake with islands" if all the boundary sites of V 
have height 2d, and if from the origin there is a path along sites having height 2d to 
the boundary. From the fact that the zeros do not percolate, and l's do percolate 
it follows that with positive probability, the origin has height 2d and is in infinitely 
many nested lakes, i.e. V\ C V2, ■ ■ ■ , V n . . ., with for i 7^ j, dVi PI dVj = 0. Consider 
a configuration from that event, consider a volume V D V n , add a grain at the origin 
and stabilize the configuration in V. In the first wave all sites will topple once because 
islands not contained in the wave would contain forbidden subconfigurations, which 
is impossible since the configuration is recurrent. After the first wave, the only sites 
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that change height are on the boundary of V n . Therefore, the origin is still unstable 
and a second wave must start. The sites included in this wave will contain the set of 
sites included in the first wave needed to stabilize r]v 7l _ 1 , but this set, with the same 
argument, is at least V n -i. Continuing like this, one sees that at least n waves are 
needed for the stabilization of rj inside V. Since with positive probability we find 
infinitely many lakes containing the origin, the number of topplings is diverging in 
the limit V \7L d with positive probability, which is what we wanted to prove. □ 

Remark 7.5. We believe but cannot prove that any translation invariant probability 
measure v G V(7i) with v(r)(0)) G (p c ,2d) is either metastable or not stabilizable. 

8 Appendix 

In this appendix we prove the ergodicity of the scenery process, used at several places 
in this, e.g. in (|4.11j) . 

Proposition 8.1. Let Q = N zd (or any space on which translations act). Suppose 
that [a is a stationary and ergodic probability measure, and X n , n G N is symmetric 
nearest neighbor random walk on Z d . Then, if initially rj is distributed according to 
Li, the process Tx n T] defined by 

T x n v(%) =v(x + X n ) 
is a stationary ergodic (in time n) Markov process. 

Proof. Let Sd denote the set of unit vectors of Z d . The process Tx n T] is clearly a 
Markov process with transition operator 



To prove ergodicity (in time n) of the Markov process, we have to show that if / is a 
bounded measurable function with Pf = f, then / is constant //-almost-surely (see 
e.g. P3.j)- So suppose that Pf = f, then, by translation invariance of li, 



Hence r e f = f, ll almost surely for all unit vectors, and hence r x f = f fi almost 
surely for all x G Z d . By ergodicity of li (under translations), this implies / = J f li 



Finally, we prove a fact about recurrent configurations that we used in our con- 
structive example. 

Define two (possibly unstable) height configurations rj, £ on V equivalent if there 
exists my : V — ► Z such that rj — £ = Ay my- It is well-known that every equiva- 
lence class of this relation contains exactly one recurrent configuration. So if rj is a 




/x-almost surely. 



□ 
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(possibly unstable) height configuration which upon stabilization yields a recurrent 
configuration £, then the relation 

j] - A v m v = f (8.2) 

simply means that stabilization of r\ requires m v (x) topplings at site x G V and yields 
£ as end result. 

Proposition 8.3. Suppose that V cZ 2 is a rectangle and rj is a recurrent configura- 
tion in V. Then upon addition of 2 grains to the corner sites and 1 grain to all other 
boundary sites, every site will topple once and the resulting configuration remains 
unaltered. 

Proof. Let 1 denote the column indexed by x G V of all ones, then (Ayl) x equals 4 
minus the number of neighbors of x in V. Therefore, the simple identity 

T] + Ay I - Ay I = T] 

shows that if rj is a recurrent configuration, addition of (Ayl) x reproduces rj, and 
makes each site topple once (cf. (|8.2jl ). 

□ 
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